Three body on-site interactions in ultracold bosonic atoms in an optical lattice and 

superlattice 
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The Mott insulator-superfluid transition for ultracold bosonic atoms in an optical lattice has been 
extensively studied in the framework of the Bose-Hubbard model with two-body on-site interactions. 
In this paper, we analyze the additional effect of the three-body on-site interactions on this phase 
transition in an optical lattice and the transitions between the various phases that arise in an optical 
superlattice. Using the mean-field theory and the density matrix renormalization group method, we 
find the phase diagrams depicting the relationships between various physical quantities in an optical 
lattice and superlattice. We also propose a possible experimental signature to observe the on-site 
three-body interactions. 
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Studies of ultracold bosonic atoms in optical lattices of- 
fer many opportunities for exploring a variety of quantum 
phases. After the first theoretical prediction of superfluid 
(SF) to Mott insulator (MI) transition in bosonic systems 
by Fisher et al. [ll , in three dimensional optical lattice by 
Jaksch et al. (2| which was followed by an experimental 
observation by Greiner et al Q, a number of quantum 
phases have been reported in the literature d, [^ . The 
ability to fine tune the control parameters in such exper- 
iments makes it possible to probe these exotic quantum 
phases [a,01- Apart from the hopping and the on-site two 
body interactions, the on-site three-body interactions can 
also influence the onset of different phases. Zhang et al. 
had earlier found the extension of the insulating lobes in 
the presence of the on-site three-body interactions, us- 
ing the decoupling mean- field theory [g. The generation 
of effective three- and higher-body interactions by two- 
body collisions of atoms confined in the lowest vibrational 
states of a three-dimensional optical lattice has been re- 
ported by Johnson et al. y\. The effect of three-body 
interactions on the insulating lobes in an optical lattice 
has been considered using the mean-field and functional 
integral approaches in the Bose-Hubbard approximation 
for optical lattices Q, UM ■ On the other hand, studies on 
optical superlattices have revealed the existence of phases 
other than the usual MI and SF phases; namely the su- 
perlattice induced Mott insulator(SLMI) phases which 
have a density modulation in the system [ll|, llJI ■ 

On the experimental side. Will et al. [Tj] have detected 
and precisely measured the on-site three and higher body 
interaction strengths by observing the collapse and re- 
vival of the superfluid matter waves in a deep opti- 
cal lattice. Nagerl et al p3| have been able to pre- 
cisely determine the on-site interaction energies including 
multi-body interaction shifts. In another work, Greiner 
et al. 15] have determined the three-body interaction 
strengths by using occupation-sensitive photon-assisted 
tunneling. However, the effect of three-body interactions 



in optical superlattices has still not been investigated to 
the best of our knowledge. 

In this paper, we study the effect of the on-site re- 
pulsive three-body interactions in addition to the on-site 
two-body interactions on various phases exhibited by ul- 
tracold bosonic atoms in an optical lattice and a super- 
lattice. We use the mean-field decoupling approximation 
and the finite size density matrix renormalization group 
(FS-DMRG) method for various densities and three-body 
interaction strengths and then compare the results. The 
system of bosons in an optical superlattice with three- 
body interaction can be described by the modified Bose- 
Hubbard model as follows:- 
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Here, aj{di) is the creation (annihilation) operator which 
creates (destroys) an atom at site i, hi — a\ai is the 
number operator, t is the hopping amplitude between 
the adjacent sites (j,j), U and W represent the on-site 
inter-atomic two-body and three-body interactions re- 
spectively, /i is the chemical potential, and A is the su- 
perlattice potential. We consider a bipartite lattice with 
sublattices A and B with a periodicity of two sites. We 
apply standard decoupling approximation [12, [la, [l3| to 
the hopping term in Eq. ([T]) to obtain the mean-field 
Hamiltonian given by 
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where the superfluid order parameter ipi —< di > is 
taken to be real [T7|, 4>i = - X)^ V'i-i-i5, the summa- 
tion over S is taken over z nearest neighboring sites. 
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FIG. 1. (Color online) Phase diagram of Eq. ([2]) for different 
W for optical lattice. The lobes represent Mott insulator 
phases for densities p — 1,2, 3, 4. 



U — U/zt, W = W/zt, JI = ^/zt and A^ = Xi/zt are 
dimensionless parameters. For an optical lattice, Ai = 
for all i, thus ipi = ip. For our optical superlattice, 
Xi = for sublattice A and A^ = A for sublattice B, 
thus V'i = 4'a{4'b) if * belongs to sublattice A {B). The 
mean-field eigen value equation is solved self-consistently 
to obtain the local superfluid density p| — V'f and density 
Pi — {fii) of the ground state of the system. 

To study the effect of W on MI phases in an optical lat- 
tice, we first present the mean- field phase diagram for an 
optical lattice (Fig. [Ij, in the U - JI plane obtained from 
the density p and the superfluid density p^, for various 
values of W. Figure [2] shows the Ji - p,p^ plot for differ- 
ent W . For an optical superlattice we show the effect of 
W on MI and SLMI phases in the phase diagram plotted 
in the fl - X plane (Figs. [3] and |4]). In Fig. [3] we present 
the phase diagram for C/ = 10 and W = 0.0. Fig. |4]is 
the phase diagram for JJ — IQ and W = 5.0. Lobes Lp 
represent the MI phase with density p. Lobes Rl to R6 
represent SLMI phases with density in sublattices A {B) 
respectively given by 1(0), 2(0), 2(1), 3(1), 3(2) and 4(2). 

The DMRG results obtained for the optical lattice are 
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FIG. 2. (Color online) Variation of density p and superfluid 
density p" with /I for [/ = 10, for optical lattice. Top to 
bottom, the first four curves represent density p and the next 
four curves represent superfluid density p". The plateaus in 
the p plots represent MI phases with vanishing p" . 



FIG. 3. (Color online) X - p phase diagram for U = 10, 
W = 0.0, for optical superlattice. 



given in Figs. [7] and [H for p = 2 and 3 respectively. 
Figures 9(a) and 9(b) are the phase diagrams for the 
optical superlattice for p — 3/2 and 2 respectively, with 
two values of W{^ 0.0, 5.0). 

For the calculations using the FS-DMRG method, we 
fix the hopping matrix element t = 1 to fix the en- 
ergy scale (so all the quantities plotted, are in units of 
t) and to estimate the critical points Uq, we perform 
finite-size scaling of the single particle gap Gl (defined 
by the difference between the energies needed to add 
an atom and remove an atom from the system). The 
plots of LGl for different system sizes L (see Fig. ([5])), 
assuming that the SF to MI transitions be long to the 
Berezinskii-Kosterlitz-Thouless (BKT) type [Illil, co- 
alesce in the superfluid phase below the Uc- The value 
of Uc is then estimated within an error bar of 0.1 if the 
values of LGl, say for L = 140 and 200, differ by less 
than 4 %. At the BKT transition the gap closes satisfy- 
ing the relation, Gl ^ exp[— a/|[/— C/c|^^^], where a is a 
constant. The correlation length ^, is finite in the gapped 
phase and diverges at the critical point as {1/Gl) — 
exp [a/\U — C/cT^^]- Near Uc, the flnite-size-scaling re- 
lation LGl[1 + {1/(2 In L + C)}] = F{^/L)^ used to 
estimate the transition point as done in Ref. [201 ■ There- 
fore, if we plot ln{L/^) vs. LGl (1 + (1/(2 In L + G)) 
then the curves for different lengths collapse in the vicin- 
ity of C/c (scaled) (see Fig. ([5]), main panel). Combining 
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FIG. 4. (Color online) X - p phase diagram for U 
W = 5.0, for optical superlattice 
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FIG. 5. (Color online) Scaling of gap LGl plotted as function 
of f/ for p = 2 and W = 0.0. 
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FIG. 7. (Color online) Phase diagram for p — 2 for various 
values of W , for optical lattice. 



the scaling method described above and DMRG results, 
we give an approximate value of Uc in various configu- 
rations. 

From the mean field results (Figs. [T] and [H) for the 
optical lattice, we find that the p = 1 MI lobe remains 
unaltered in the presence of W . However, for higher den- 
sities, the critical value Uc{W) for SF-MI transition de- 
creases as W increases (e.g. Uc for p = 2 lobe decreases 
from ^ 10 to ^ 7 when W increases from 0.0 to 4.0) and 
this is more prominent as the density increases, as shown 
in Fig. [TJ Also, the MI lobes get enlarged as the W in- 
creases. This is further confirmed by the DMRG results 
from Figs. [7] and [5] The trend is the same for both the 
densities, p — 2 and 3, but the effect of W is more when 
the density is large. For p — 2, UciW — 4.0) ~ 1.6 com- 
pared to UciyV — 0.0) « 5.7. For p = 3, Uc decreases 
steadily as W increases; Uc{W = 0.0) « 8.6, Uc{W = 
1.0) « 6.6, UciW = 2.0) « 4.6, for UciW = 3.0) « 2.6 
and Uc{W = 4.0) w 0.8. The reason for this behavior 
at higher densities is that there is a greater probability 
of having three or more atoms at a site, which enhances 
the three-body interaction and suppresses atom hopping 
from one site to another. 
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FIG. 6. (Color online) Scaled gap Gl plotted as a function 
of scaled U for p = 2 and W = 0.0. The curves for different 
lengths collapse in the vicinity of Uc as the correlation length 
^ diverges exponentially near Uc- (Inset) Scaled gap Gl plot- 
ted as a function of U. The curves for different lengths cross 
at Uc (~ 5.8) showing the critical point for SF-MI transition. 



From the mean field results for the optical superlattice 
(Figs. [3] and H]), we see that the lobes LI, Rl and R2 
remain unaffected in the presence of a finite W. This is 
expected because in such configurations, no two adjacent 
sites have more than two atoms, and for an atom to hop 
two sites is a second order process, which is of much less 
probability. However, the SLMI phase R3, which has 
sublattice atomic densities pA = 2 and ps = 1 (and 
thus has average density p — 3/2) gets enlarged in the 
presence of W. This is understood from the following 
reason. When W = 0.0 and as we increase A, keeping 
the average density p = 3/2, the ground state goes from 
the superfluid to the SLMI phase R3. With a further 
increase of A > 15, the ground state is again a superfluid 
which has 2 < p^ < 3 and < ps < 1. However, 
when W is finite, the system does not prefer to have 
sublattice densities above 2. Thus SLMI phase R3 has a 
lower energy than superfiuid with sublattice densities 2 < 
PA < 3 and < PS < 1. DMRG resuhs [Fig.9(a)] show 
a similar trend. When W = 0.0, the system undergoes 
a phase transition from the SF phase to the SLMI phase 
R3 at a value of A ~ 0.3. However, in the presence of 
a finite W, the transition occurs at a lower A(^ 0.15), 
signifying the enlargement of the insulating lobes. 

The phase diagram for the system with a filling factor 
p = 2, shows a marked difference in the presence of W. 
Comparing Figs. [3] and HI we find that the MI lobe L2 




FIG. 8. (Color online) Phase diagram for p = 3 for various 
values of W, for optical lattice. 
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FIG. 9. (Color online) (a) Phase diagram for p — 3/2 and 
(b) phase diagram for p = 2, in an optical superlattice with 
W = 0.0 and 5.0. 



becomes large and its tip shifts from A ^ 1 to A ^ 9. 
Also as A increases to -- 17.5 (Fig. ^, the MI lobe L2 
goes to the SLMI phase R4. However, in the presence of 
W (Fig. lU, the tip of the R4 lobe gets shifted to ~ 20.5. 
As we have considered the maximum value of A till 18.0, 
SLMI phase R4 does not appear in Fig. S) Similar results 
are also obtained by the DMRG analysis (Fig. 9(b)). For 
finite W = 5.0, the critical superlattice potential Ac for 
transition from the MI phase L2 to SF phase shifts from 
9.4 to 14.8 and that for the SF phase to SLMI phase R4 
shifts from 10.6 to 15.2. In the absence of W and for 
lower values of A, the superlattice initially stays in the 
MI phase L2, for U = 10.0. As A becomes comparable to 
U, the system goes from the MI phase to the SF phase at 
A ~ 9.4. As A is further increased, the system goes from 
the SF phase to the SLMI phase R4, at a value of A ^ 
10.6. For W = 5.0, the system initially is in the MI phase, 
L2. Now, due to the presence of W, the SF window is 
shifted to a A value which is comparable to U + W a,s 
shown in Fig. 9(b). The MI to SF transition takes place 
at A ^ 14.8, and the second transition from SF to SLMI 
(R4) takes place at A ~ 15.2. The SF window not only 
shifts for M^ = 5.0 but also shrinks when compared to 
that oi W ^ 0.0. This shifting of the R4 lobe can be 
understood as follows: in the R4 phase there are three 
atoms at every alternate site. As W increases, it becomes 
difficult to confine three or more atoms at a single site. 
Hence to suppress this effect due to the increase in W, 
we need to increase A. 

The three-body interaction strength scales with the 
two-body interaction strength as follows: W oc 
ln(C7y2)(l/o/£;,)3/4 e-^V^°^^alPU^ i MM- The 
typical range of a^fc^ is 10^* to 10^^ which supports 



the fact that the three-body interaction is weaker than 



the two-body interaction [2j|. Also it can be seen that 



three-body interaction is tunable and can be adjusted 
by varying (Vo/Er). The three-body effects have been 
experimentally observed before through various methods 
as mentioned earlier [13l - [l5| . We propose an alternate 



method to observe these effects in an optical lattice and 
superlattice that we have considered in our present work. 
The effect of W is very small compared to the two-body 
interaction in the system of bosons in an optical lattice. 
This is of course true when the filling factor of the sys- 
tem is unity. From Eq( [T]) it is clear that the three-body 
energy scales as n^. Therefore, in order to observe the ef- 
fect of the three-body interaction in the experiment it is 
important to study the SF-MI transition at higher densi- 
ties. In the seminal work of Greiner et al. Q, the SF-MI 
transition was observed by probing the excitation spec- 
trum resulting from a particle-hole excitation. Such an 
excitation was created by applying a potential gradient 
to the system in the MI phase. By plotting the excitation 
probability versus an applied vertical potential gradient, 
two narrow resonance peaks were seen. The first peak 
was at the potential gradient equal to the single particle 
excitation gap, and this corresponds to the MI shell at 
density equal to one. One of the possible reasons for the 
appearance of the second peak was the particle-hole ex- 
citation created in the MI shell at a density equal to two. 
In the MI shell at a density equal to two, the particle- 
hole excitation at a given site would populate one of the 
neighboring sites with three-atoms. In principle, when 
there are three or more atoms in a lattice site, the atoms 
will experience the effect of W along with that oi U. In 
general when there are n (> 1) atoms in each site the 
system is in the MI phase with a density equal to n, the 
excitation gap is A = U + {n—l)W for the optical lattice 
and A = U + {n — 1)W + A for the optical superlattice. 
Therefore, by measuring the values of the potential gradi- 
ents for the higher order peaks, and taking the difference 
between them for different densities, it would be possible 
to determine the value of W. 

In conclusion, we have studied the effects of the onsite 
three-body interactions in a system of neutral bosons in 
an optical lattice and superlattice. We first use the mean 
field theory to understand the behavior of the system 
and then confirm the results using the DMRG method. 
In the optical lattice, and the superlattice as well, we find 
that the Mott insulator lobes get enlarged as the value 
of W increases. When the density p = I, the effect of 
W is not significant. However, as the density of the sys- 
tem increases the effect of W becomes significantly large 
which changes the SF-MI critical point drastically. We 
obtain the phase diagrams for different combinations of 
densities, strengths of the three-body interaction and the 
superlattice potential. Finally, we have also suggested a 
possible experimental scenario by which it may be possi- 
ble to observe a signature of the three-body interaction. 
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